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ARRANGEMENT OF FUSIFORM BODIES TO REDUCE THE 
WAVE DRAG AT SUPERSONIC SFEEDS 1 
By Morris D. Friedman and Doris Cohen 


SUMMARY 


Using linearized slender-body theory and reverse -flow theorems, the 
wave drag of a system of fusiform bodies at zero angle of attack and 
supersonic speeds is studied to determine the effect of varying the rela- 
tive location of the component parts. The investigation is limited to 
two-body and three-body arrangements of Sears-Haack minimum-drag bodies. 
It is found that in certain arrangements the interference effects are 
beneficial, and may even result in the two- or three-body system having 
no more wave drag than that of the principal body alone. The most favor- 
able location appears to be one in which the maximum cross section of 
the auxiliary body is slightly forward of the Mach cone from the tail of 
the main body. The least favorable is the region between the Mach cone 
from the nose and the forecone from the tail of the main body. 


INTRODUCTION 


When an airplane is to be equipped with external fuel tanks or 
prominent nacelles, the effect on the drag will vary widely with the 
location of such auxiliary bodies relative to the other parts of the 
airplane. In reference 1, calculations were made of the theoretical 
interference drag between the fusiform bodies of some typical arrange- 
ments, under the conditions of supersonic speed and zero angle of attack. 
Later developments in linear theory have provided a simpler method of 
performing such calculations, and the present paper is a revision of 
reference 1 to take advantage of these developments. Both reference 1 
and the present work are largely based on suggestions of R. T. Jones. 

Two arrangements will be considered - a two -body combination, as 
when one body is suspended beneath another, and a laterally symmetric 


■'■Supersedes NACA RM A 51120, "Arrangement of Bodies of Revolution in 
Supersonic Flow to Reduce Wave Drag" by Morris D. Friedman, 1951 « 
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three -body arrangement. The radial and streamwise displacements of the 
auxiliary body or bodies relative to the main one will constitute the 
parameters of the investigation. The calculations will be made for com- 
binations of Sears -Haack minimum-drag bodies (refs. 2 and 3), but the 
method of analysis is applicable to any slender shapes for which the 
pressure fields are known. In particular, it may be mentioned that the 
main body and auxiliary bodies need not be similar. 


ANALYSIS 

Reversed-Flow Theory 


A basic condition of the analysis is that the bodies be slender 
enough so that they may be represented by a linear distribution of sin- 
gularities - sources for a body of revolution, or higher-order singulari- 
ties for cambered bodies - of which the strength may be determined from 
local conditions. In that case, simultaneously reversing the direction 
of flow and the sign of the source strength associated with a given iso- 
lated body does not change the shape of the body, and the reversed-f low 
theorems of reference 4, which are stated in terms of source distribu- 
tions, may be applied to the bodies themselves. However, the stream-lines 
at a distance from the body are altered, so that it is not generally to 
be expected that the theorems would be applicable to a system of bodies 
of prescribed geometry. The location of the individual bodies of the 
system in the streamlines of the other bodies has the effect of introduc- 
ing additional camber into the boundary conditions and thereby modifying 
the equivalent distribution of singularities. Calculations made to inves- 
tigate the effect of such induced camber on the drag of slender bodies 
indicate that the magnitude is not likely to be any significant fraction 
of the thickness drag. 2 The additional drag introduced by the induced 
camber will therefore be ignored. With this simplification, the drag of 
a system of slender bodies may be said to remain unchanged when the 
direction of motion is reversed. 

In the present calculations, in which each of the bodies is symmet- 
rical fore and aft, the first consequence of the reversibility property 
is that only rearward (or forward) displacements of the auxiliary bodies 
relative to the main body need be considered. 

The reversibility property also leads to the possibility of com- 
bining the pressure fields for forward and reverse flow before computing 
the drag, and taking half the result as the drag in either direction. 


2 For example, the drag due to parabolic camber of a line of sources 

corresponding to a Sears-Haack body is ^ (M 2 - l) h 2 times the drag of 

of the uncambered body, h being the maximum camber in percent of body 
length. 
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This point may he demonstrated as follows: 

Let dR/dx he the local inclination of any element of surface area 
(of a body of revolution) to the stream and Cp^ he the pressure coef- 
ficient at the centroid of the element when the body is in forward motion. 
Then the corresponding element of drag (sketch (a.)) is simply (to the 
first order) 

q cp^. dR/dx dS 

where q is the free -stream dynamic 
pressure and S is the surface area. 

Now, let the body he reversed on the 
x axis (sketch (b)) and flown tail 
foremost. At the previously con- 
sidered element of area a new pres- 
sure coefficient Cp^ will result; 
the slope dR/dx will merely he 
reversed in sign. The sum of the 
two elements of drag will he 

+ “r - q(c Pf - =p r ) (n) f 48 

( 1 ) 

where the subscript f refers to _ 

quantities measured in forward V D r 

motion and r to quantities in 
rearward motion. The total com- 

hined drag is the integral of this Sketch (h) 

quantity, and is twice the drag of the body traveling in either direction. 

The foregoing device, similar to one first suggested in refer- 
ence 5, results in considerable mathematical simplification in many 
cases. Thus, the argument of reference 6 that in the combined flow 
field the drag of one wing due to the field of another is equal to that 
of the second due to the field of the first can be extended to apply to 
systems of slender bodies under the conditions outlined above. Then 
only one calculation of the interference drag need be made for each pair 
of bodies; for example, in the present analysis only the drag of the 
auxiliary body due to the combined pressure field of the main body will 
be calculated. 
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Pressure Field of a Sears -Haack Body- 


In the examples to he worked in this paper, all bodies will be 
closed bodies of revolution having individually the form for minimum 
theoretical wave drag for given length and volume. 3 This shape (refs. 2 
and 3) is given by 


R 2 (x) = c(Z 2 - x 2 ) 3 ^ 2 


( 2 ) 


(See Appendix for a list of symbols.) The flow about one such body can 
be calculated as the flow due to a distribution of sources and sinks 
along the body axis, the source strength f (x) being related to the body 
geometry by the equation (ref. 7) 


2 it f ( x) 




( 3 ) 


The pressure coefficient near the body is given to the first order 
(see ref. 8) by the relation 


C P - - Y -(f) • (M 

in which u = ckp/dx is the streamwise component of the perturbation 
velocity and v = c)cp/c)r the radial component. Inasmuch as the latter 
component falls off with distance from the body as l/r and its con- 
tribution to the pressure decreases with l/r 2 , it may be possible to 
neglect the second term in Cp in computing the interference drag. 

It will be shown by a numerical example that this simplification is in 
fact permissible in the present investigation. At the surface of the 
body, the ratio v/v is, to first order, the streamwise slope of the 
body and, in the. case of the symmetrical bodies under consideration, 
the integration for the drag will result in canceling out all effects 


3 In a recent analysis, Conrad Renneman, Jr., of the NACA has found 
that the form of a body of revolution for minimum drag in the field of a 
larger body is not significantly different from the Sears -Haack shape, 
nor can any important reduction in the theoretical wave drag be achieved 
by modifying the shape of the body from the Sears -Haack shape. 
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of the radial component of velocity. Thus it is sufficient for our 
purpose to retain only the linear term in equation (ha .) , writing 



(>«>) 


From reference 7, 


cXP = _ r 
<3x J 


x-pr f'(£)d£ 


nose 


J (x-£) 2 - P 2 r 2 . 


(5) 


From equations (2) and (3), f'(l) is in the present case 


3 Z 2 -2| 2 

^ Vc 


<\/z 2 - l 


2 


and 


Cr> ~ -3c 


/ 


*>3 r (l 2 - 2| 2 )d£ 


_ z y ( Z 2 -! 2 ) [(x-|) 2 - P 2 r 2 1 


( 6 ) 


It should be remarked that for £ > Z the integrand vanishes, and that 
Cp is zero when x - Pr < - Z. 

The integration in equation (6) yields three different expressions, 
depending on whether 

pr-Z<x<Z-pr (Region I) . 

x - Pr < .Z < x + Pr (Region II) 


or 


Z < x - Pr 


(Region III) 


The three regions defined above are shown in figure 1. The expressions 
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for the approximate pressure coefficient are ' 


Region I 



Z 2 +2(3r( Z-x) -2x 2 
J (Z+pr) 2 -x 2 


K 0 (ki)-<\/ (Z+Pr) 2 -x 2 Eo(k!) + 2xA c (t 1 ,k 1 ) 


(Z-pr) 2 -x 2 
( Z+pr) 2 -x 2 


kn = 


♦i 


sin -1 



Z+pr+x 

21 


(7) 


Region II 


c 


II 




-3«c 


Z( Z+2pr-2x) 


2 



K Q (k 2 ) - 2/z^ r E Q (k 2 ) + 2xA 0 (* 2 ,k 2 ) 


k 2 - 


x 2 -( (3r-Z) ; 
4zpr 


\|f 2 = sin 


-/ 


2Z 


Z+pr+x 


(8) 


Region III 


c p ~ 3nc 
III 


(x-pr) ; 


J x 2 -(pr-Z) : 


K Q (k 3 ) + J x 2 -(pr-Z ) 2 E 0 (k 3 ) - 2xA 0 (t 3 ,k 3 ) 


k 3 



4Zpr 


■Or-Z) ; 


*3 


sin -1 



Z+x-Pr 
Z+x+Pr , 


(9) 


The quantities K Q , E Q and A 0 ' are defined in Appendix A and tabulated 
in reference 9 ; A Q is tabulated also in reference 10. 

The magnitude of the approximate pressure coefficient at the sur- 
face of a body of fineness ratio 10, as given by equation (7), is plotted 
in figure 2(a) , together with the more accurate values obtained by the 
use of equation (4a). It is apparent that the difference, in the case 
of so slender a body, is not great and, as previously noted, will diminish 
rapidly with distance. 

An isometric sketch of the pressure coefficient calculated by equa- 
tion (4b) is shown in figure 2(b). Of particular interest is the 
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logarithmic infinity along the Mach cone from the tail of the body. 
Except at the body itself, the pressure is finite everywhere else and 
goes smoothly through the forecone from the tail, in spite of the change 
in the form of its mathematical expression. 


Combined Pressure Field 


If now the body is reversed in heading and the resulting pressure 
field subtracted from that given above, the combined pressure field is 
found as 

c p (x) = c p (x) - Cp(-x) (10) 

because of the fore-and-aft symmetry of the body. The various regions 
of the combined field are shown in figure 3* 

It is seen from figure 3 that there is only one region in which any 
possiblity of further mathematical simplification appears. In this 
region (where Region I and Region I of the reversed field overlap) we 
have to consider 

Cp = c p ( x ) ” c p ( -x ) (ll) 

I P I *1 

Through the relation (ref. 10, p. 36) 


when 


A 0 (i|r,k) + A 0 (\jf,k) = 1 + k ,2 K Q sint sint 


tan tan i|r 


1 

k 


equation (ll) reduces to 



-6«cx 


( 12 ) 


Thus, the pressure gradient of the combined field is a constant in the 
neighborhood of the body, as specified in reference 6 for minimum wave 
drag with a given volume. 4 

4 A1 though reference 6 deals specifically with thin wings, it is 
readily shown that the same considerations hold for slender bodies under 
the assumptions used herein. 
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Calculation of the Wave Drag 

As previously indicated, it is proposed in calculating the drag to 
ignore the drag introduced by the curvature of the flow due to adjacent 
bodies, and therefore to replace each body by the equivalent source 
distribution in a uniform stream. If, following this course, we con- 
sider the entire flow field to be essentially the result of the super- 
position of the fields of the individual bodies, we have to compute 

(a) one-half the drag of each body in its own combined pressure 
field, 

(b) the drag of each auxiliary body in the combined field of the 
main body, and 

(c) the drag of one auxiliary body in the field of the other, if 
more than one auxiliary body is included. It may be seen that the drag 
of the second body in the field of the first is taken care of by the 
factor of 2 introduced by the use of the combined pressure field. 

If equation (l) is applied to a body of revolution, the total drag 
in combined flow may be written 


dx 


and, for the Sears-Haaek bodies. 


■N+ 1 


D = 2:rtq 


Z 


? R dR 
C P R dx 


= 


- dB 2 

P ^T 


dx 


(13) 


-Z 


dR 2 

dx 


= -3cx J l 2 -x 2 


(14) 


The drag of each body due to its own pressure field is then one -half 
that obtained by substituting from equations (12) and (l4) in equa- 
tion (13) , or 


Do = 9* 2 qc 2 


^ + 1 


x 2 Z 2 -x 2 dx 


-Z 

= | jt 3 qZ 4 c 2 

8 


( 15 ) 


or, since the maximum cross-sectional area is, from equation (2), 


Sp> = jtcZ 3 


( 16 ) 
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the drag in terms of frontal area is 



which is in agreement with the value given in reference 2. 


(17) 


The interference drag is calculated by substituting in equation (13) 
the pressure coefficient associated with one body and the value of dK 2 /dx 
associated with the other. Considering first a two-body combination, we 
may take the origin of coordinates at the center of one, which will be 
designated the "main" body, and let the center of the auxiliary body be 
displaced from it a distance x Q downstream and a distance r Q laterally 
or vertically. Then, if for a first approximation the pressure field of 
the main body is assumed not to vary significantly in the distance between 
the axis of the auxiliary body and its surface-, the values of c p are 
obtained from equations (12), (8), and (9) by letting r = r Q . Equa- 
tion (l4) is modified to take into account the displacement of the auxil- 
iary body: 

2 

-3c 1 (x - x Q ) - (x - x 0 )^ (18) 

and the geometric characteristics of the second body are used to deter- 
mine c x and l x • 

It is immediately apparent that if the auxiliary body lies entirely 
within Region I (fig. 3), the value of c p , and therefore of the drag, is 
independent of r Q . In fact, since the pressure gradient is a constant 
throughout the region, the drag is entirely independent of position so 
long as the body remains within Region I. The interference drag in this 
case is simply 


x D 1 = | * 3 (19) 

or, more generally, 'D 1 /q equals the pressure gradient 6itc times the 
volume of the auxiliary body, regardless of its shape. 

In Region II, however, the pressure gradient changes from negative 
to positive (see fig. 2(b)) and a small body placed so as to take advan- 
tage of this buoyancy would conceivably experience a negative interfer- 
ence drag, or thrust, which would act to reduce the total drag of the 
combination. Substitution of c p and, when required, c p in equa- 
tion (13) results in integrals which can be evaluated only numerically. 
We therefore proceed at this point to the. consideration of numerical 
examples . 
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NUMERICAL CALCULATIONS 


For an exploratory investigation, the simple case of two bodies of 
fineness ratio 10, the small body having one -half the length of the main 
body, was chosen. This combination is derived by assigning the following 
values to the parameters of the problem: 

c = 0.01 l = 1 

c x = .02 l x = l/2 

In compressible flow, the Mach number and cross-stream dimensions always 
enter together in the form of 0r. Three values of 0r Q /l have been 
selected: 0.25, 0.5* and *1.0. The effect of streamwise displacement is 

investigated to a distance equal to the length of the larger body. 

The interference between the two outer bodies of a three -body con- - 
figuration was calculated in the case of 0r o /l = 0.25, the only case, 
because of the limited zones of influence at supersonic speeds, in which 
any such interference takes place. The interference drag was found to 
be negligible. The results to be presented are therefore equally appli- 
cable to the two -body or three -body arrangement, only a factor of two in 
the interference drag being required. 


RESULTS AND DISCUSSION 


The variation of the interference drag with streamwise and radial 
displacement is shown in figure 4. Because of the symmetry of the curves, 
the effects of rearward displacement only are presented. The sketches 
indicate, for M = V~2, the relative positions of the bodies at which the 
interference is greatest. The anticipated favorable interference is 
observed when the small body is situated astride the region of negative 
pressure gradient just ahead of the Mach wave from the stern of the large 
body. A second case of favorable interference is seen in figure 4(c) 
when the pressure field of the large body acts only on the rear of the 
small body, resulting in unopposed thrust. This situation, in the case 
of forward displacement of the small body, is also responsible for the 
dip in the interference drag at x Q /l = 0.7 when 0r Q /l =0.5, and rein- 
forces the minimum in the neighborhood of x 0 /l = 0.95 when 0r o /z = 0.25. 5 


The situation described actually takes place at x 0 /l = -0.5 
and -0.75 in the two cases, but the general slope of the curve due to 
the proximity of stronger minimums causes a shift of the secondary dips. 
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In the corresponding cases of rearward displacement, the thrust is the 
result of the location of the rear of the large body in a positive 
gradient due to the small body, but there is a more involved balance of 
forces and the net result is not as easily foreseen. As is to be 
expected, the maximum benefit obtainable decreases with increasing radial 
separation of the bodies. The maximum penalty is incurred when the 
auxiliary body is added between the Mach cone from the nose and the fore- 
cone from the tail of the main body. 

In figures 5 and 6 the wave drag coefficient based on total frontal 
area is shown for two-body and three-body arrangements. These figures 
indicate that in frictionless potential flow it might be possible to 
increase the volume by as much as 25 percent (three-body arrangement) 
and at the same time actually decrease the wave drag. In practice, of 
course, the additional friction drag might easily nullify any such gain. 
Nevertheless, if there are to be auxiliary bodies, the importance of 
their relative positions seems clear. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Sept. 8, 1954 
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APPENDIX 

SYMBOLS 

p2 

n ■ 

coefficient containing "body dimensions, — g — 

Z 

value of c for auxiliary body • 

% 

wave-drag coefficient, -~- 

- qSp 

Di ■ 

interference drag coefficient, 

qSp 

local pressure coefficient 

pressure coefficient in forward motion 

pressure coefficient in reversed motion 

pressure coefficient in "combined" motion, Cp_^ - Cp_^ 
wave drag 

wave drag of isolated body- 
interference drag 
wave drag in forward motion 
wave drag in reversed motion 
combined wave drag, Df + D r 
2 

— times complete elliptic integral of the second kind 
source strength 

modulus of elliptic integrals (with subscripts to indicate 
different values) 

complementary modulus, \/~l ~ k^ 

2 

— times complete elliptic integral of the first kind 
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l 

h 

M 

q 

r 

r o 

R 

^max 

S 

S c 

&p 

u 

V 

V ' 

X 

x o 

p 

A 0 (t,k) 

t 

9 

♦ 


half-length of body 

half-length of auxiliary body 

free -stream Mach number 

free-:stream dynamic pressure 

radial coordinate, measured from body axis 

radial coordinate of center of auxiliary body 

local radius of body (function of x) 

radius of body at maximum cross section 

surface area 

cross-sectional area 

maximum cross section, or frontal area 
streamvri.se component of perturbation velocity 
radial component of perturbation velocity 
free -stream velocity 
streamwise coordinate 

streamvri.se coordinate of, center of auxiliary body 

-v/m 2 - 1 

Heumann's elliptic function, E 0 (k)F(ty,k*) + K 0 (k)E(t,k l j - 
K 0 (k)F(t,k' ) , tabulated in references 9 and 10 

streamwise coordinate of source 

.perturbation velocity potential 

argument of incomplete elliptic integrals in A 0 (with sub- 
scripts to denote various values) 
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Figure /.- Regions in pressure field around isolated Sears - Haack body. 






on body 





18 


NACA TN 33^5 





4 


Figure 3.- Various regions in combined pressure field of a single body . 










Figure 5. - Wave - drag coefficient of two -body systems based on total frontal 








